arXiv: 1509.06463v2 [math.NA] 25 Feb 2016 


MATHEMATICS OF COMPUTATION 
Volume 00, Number 0, Pages 000-000 
S 0025-5718(XX)0000-0 


FINITE ELEMENT EXTERIOR CALCULUS WITH 
LOWER-ORDER TERMS 

DOUGLAS N. ARNOLD AND LIZAO LI 


Abstract. The scalar and vector Laplacians are basic operators in physics 
and engineering. In applications, they frequently show up perturbed by lower- 
order terms. The effect of such perturbations on mixed finite element methods 
in the scalar case is well-understood, but that in the vector case is not. In this 
paper, we first show that, surprisingly, for certain elements there is degrada¬ 
tion of the convergence rates with certain lower-order terms even when both 
the solution and the data are smooth. We then give a systematic analysis of 
lower-order terms in mixed methods by extending the Finite Element Exterior 
Calculus (FEEC) framework, which contains the scalar, vector Laplacian, and 
many other elliptic operators as special cases. We prove that stable mixed 
discretization remains stable with lower-order terms for sufficiently fine dis¬ 
cretization. Moreover, we derive sharp improved error estimates for each indi¬ 
vidual variable. In particular, this yields new results for the vector Laplacian 
problem which are useful in applications such as electromagnetism and acous¬ 
tics modeling. Further our results imply many previous results for the scalar 
problem and thus unifies them all under the FEEC framework. 


1. Introduction 

The vector Laplace equation, and, more generally, the Hodge Laplace equation 
associated to a complex, arise in many applications. The discretization of such equa¬ 
tions is a basic motivation of the Finite Element Exterior Calculus (FEEC) [Ull]. In 
many applications the equations include variable coefficients and lower-order terms. 
While the former is included in the standard FEEC framework through weighted 
inner products, the latter is not, which is the subject of this work. One might 
expect that lower-order perturbations degrade neither the stability nor the conver¬ 
gence rates of stable Galerkin methods. However, this need not to be true. While 
stable choices of finite elements for the unperturbed Hodge Laplacian remain stable 
for the perturbed equation, we find that certain lower order perturbations result in 
decreased rates of convergence. Other choices of element pairs or perturbations do 
not lower the convergence rate. The situation is subtle. 

First, to fix ideas, we consider a simple example taken from magnetohydrody¬ 
namics m Chapter 3]: given vector fields /, n on a domain O C , find a vector 
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field B satisfying: 

curl curls — curl(i; x B) = f, div B = 0, in 
B ■ n = 0, (curls — vxB)xn = 0, on dfl. 

Physically, S is the non-dimensionalized magnetic field inside a conductor moving 
with a velocity field v. The system admits a solution only when div/ = 0. In 
that case, the solution also satisfies the following vector Laplace equation, which is 
solvable for any data and hence more suitable for discretization: 

(1.1) — grad div S + curl curl S — curl(i; x S) = /. 

For a mixed method, we introduce a = cuvlB — v x S and solve the coupled system: 

a — curl S + n X S = 0, curl cr — grad div B = f. 

A common stable choice of mixed elements, at least when ?; = 0, seeks S in the 
space of Nedelec face elements of the second kind of degree r > 1, and a in the 
space of Nedelec edge elements of the second kind of degree (r + 1) [T31 [T3]. In this 
case, for the unperturbed problem, that is when v = 0, the convergence for the L^- 
error in a is of optimal order if the solution is smooth enough. However, 

as we show in Section |3] and verify by numerical computation in Section [3 when v 
does not vanish, the L^-convergence for a is reduced to order 0(h’'+^). A similar 
phenomenon was observed in mixed methods for the scalar Laplacian by Demlow [B] . 
But the vector case we study has more surprises. For example, consider the vector 
Laplacian perturbed by a zeroth-order term: for some real coefficient A, 

— grad div u + curl curl u + Au = /, in H, 
u X n = 0, div m = 0, on 9H. 

This problem arises, for example in electromagnetism where u is the electric field 
and A is the conductivity coefficient. If we use the same mixed finite element 
method just considered, here with a = curlu, then the error in a is one order 
suboptimal for a general matrix coefficient A, but optimal if A is a scalar coefficient. 

We now summarize convergence rates derived from our main abstract theorems 
applied to the perturbed Hodge Laplace problem, of which the previous vector 
Laplace problem is an instance. These are important in and directly relevant to 
many applications. In particular, our results for the vector cases, that is, 1-forms 
in 2D and 1- and 2-forms in 3D, are new. On a domain H in R", a fc-form has 
(^) coefficients and the meaning of the exterior derivative d and the codifferential 
S depends on k. For example, when n = 3 and fc = 2, we have the case considered 
before with Sd + dS = — grad div + curl curl. For any 0 < fc < n, the fc-form Hodge 
Laplace equation seeks a fc-form u satisfying: 

(1.2) Lqu := {dS + 6d)u = /, 

with proper boundary conditions. Our abstract theory applies to equation (1121) 
perturbed by general lower order terms: 

Lu = [(d -|- ^ 3 )(d -|- I 2 ) -l- (d -l- l 4 ){d + li) + l^]u = /, 

with proper boundary conditions. In Section [3l we allow U to be general linear 
operators, but here we assume that they are multiplication by smooth coefficient 
fields. For example, h takes a /c-form to (fc -I- l)-form, thus may be viewed as a 
(fc^i) (fc) matrix field. 


FEEC WITH LOWER-ORDER TERMS 


3 


The mixed formulation solves simultaneously for u and for a = (6 + l 2 )u, which 
is a (fc — l)-form. On a simplicial triangulation of O, for any r > 1, we have four 
canonical pairs of mixed finite elements for (cr, u): 

Vr+iA'^~^ X X VrA'^, VrA^-^ X iP"A'=, iP-A'=-i X V-A^. 

In dimension < 3, all these are classical mixed finite elements. For example, the pair 
consisting of Nedelec face elements of the second kind and Nedelec edge elements 
of the second kind before is Vr+iA^ x VrA^. For more, the Periodic Table of 
the Finite Elements (http://femtable.org/) collects all these elements and their 
correspondence to classical elements. A main result in FEEC is that the above 
four pairs lead to stable mixed finite element methods for the unperturbed Hodge 
Laplace problem. Eurther, the rates of convergence for the L^-errors in cr, dcr, u, 
and du are optimal, as determined by the approximation properties of the spaces. 
Thus, for example, if the pair Vr+i A'^-i X VrA^ is used, the order of convergence 
for the L^-error in cr, da, u, and du are r + 2 , r + 1 , r + 1 , and r, respectively, while, 
for V-A^-^ X V-A^ all four -errors converge with order r. 

Eor the perturbed system, our discrete stability result fTheorem l3.2p implies that 
when the perturbed problem is uniquely solvable at the continuous level (which is 
the generic case, as we will prove), then the mixed discretization of this problem 
using any of the four stable pairs before is still stable for sufficiently hne discretiza¬ 
tion. Eurther, if we have full elliptic regularity (for example on a smooth domain 
m), then our improved error estimate (Theorem 13.411 implies that the con¬ 
vergence rates for the unperturbed Hodge Laplacian still hold for the perturbed 
problem with a few exceptions, as summarized in Table [I] 

Table 1. L^-error rates for EEEC elements solving Hodge Laplace equation 


elements 

a 

da 

u 

du 

Vr+lA^-^ X VrA^ 

jr + 2 , if no I 2 , k, k, 

j 

r + 1 , if no k, 

r -I- 1 

r 


lr-|-l, otherwise. 


r, otherwise. 



V-+^A^-^ X VrA^ 

r -1- 1 

i 

r -f 1, if no ^ 4 , 

r -I- 1 

r 




r, otherwise. 



VrA’^-^ X V-A^ 

J r -1- 1, if no I 2 , k, 

r 

r 

r 


1 r, otherwise. 




V-A^-^ X V-A^ 

r 

r 

r 

r 


Eor example, in our first example, H is a 2-form and equation dm) has an 1 2 
lower-order term, which leads to reduced L^-error rates for a. In practice, this 
suggests that the smaller space should be used for a = curlH, because the 

use of the bigger space Vr+iA^ does not improve the convergence rates for cr or any 
other quantity. In our second example, the zeroth-order term is a generic I 5 term, 
hence it also degrades the L^-error rate in cr. We observe that some lower-order 
terms, namely li and I 3 terms, do not degrade the error rates in any of the cases 
above and that the L^-convergence rates for u and du are unaffected by the lower- 
order terms. We also note that, most surprisingly, the lower-order term Z 4 has the 
worst effect on the convergence rates of the error in cr, degrading the L^-error rates 
in both a and da, yet a has no apparent dependence on I 4 . 
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Historically, the effect of lower-order terms on the convergence of finite element 
methods was first studied by Schatz [18], for the scalar Laplace equation with 
Lagrange elements. The key tool in that analysis was the Aubin-Nitsche duality 
argument, which was introduced to prove L^-error estimates in both non-mixed 
m and mixed methods |9]. These ideas guide the current techniques. Directly 
relevant to this work are the studies by Douglas and Roberts [Hi and Demlow 
on mixed finite element discretization of the scalar Laplace equation, which 
is the Hodge Laplace equation for n-forms in n dimensions. The primary mixed 
finite elements for this problem are the Raviart-Thomas (RT) family [T7j 

and the BDM family P^A""^ (5] |4]. Douglas and Roberts proved the optimal 
convergence rates for both variables when the RT family is used. Demlow showed 
that for BDM elements, even for constant coefficients and smooth solutions, there is 
degradation of the convergence rate for a in the problem — div(gradu-|-6u)-|-CM = / 
while there is no degradation if the same problem is formulated as — div grad u + b- 
gradu + cu = f. All these classical results on non-mixed and mixed finite elements 
for scalar Laplace problem with lower-order terms can be read off directly from 
Table [T] Our approach here gives a uniform derivation of all classical L^-estimates 
for the scalar/vector Laplacian perturbed by lower-order terms under the extended 
abstract FEEC framework. 

The rest of this paper is organized as follows. We first briefly review the basic 
FEEC framework for the unperturbed abstract Hodge Laplacian in Section[5| Then 
we lay out our extended abstract framework and state our two main discrete result: 
stability theorem (Theorem 13.21) and improved error estimates (Theorem 13.41) in 
Section [3l After that, we prove the well-posedness theorems at the continuous 
level in Section |4| Then we prove the two main discrete results in Section [5| and 
Section [5] respectively. Einally, in Section 0 we show that the estimates are sharp 
for the Hodge Laplacian case through numerical examples. 

2. Review of the abstract FEEC fraemwork 

FEEC is an abstract framework for analyzing mixed finite element methods m- 
A Hilbert complex (kF^, is a sequence of Hilbert spaces and closed densly- 
defined linear operators —>■ with closed range satisfying d^+^od^ = 0. 

We use (•, •) to denote the kF-inner product and || • || to denote the kF-norm. Let 
dfc_|_i be the adjoint of d*, = D(d^), and := From this definition, 

{d^u,v) = (u, d^_|_;^z;), for all u £ and v € Ffe+i- One important structure is the 
Hodge decomposition: 

kF'^ = 3'= © (3'=)-l^ = »'' © ^'^ © (3'=)-^'^ = »'= 0 

where 5'^ = kerd'^, = imd^-\ = 3^= n and S* = imd*+i = (3'')-^^. 

We assume the Hilbert complex satisfies the compactness property, that D 
is a compact densely embedded subspace of kF. This is true for all the cases we 
are interested in. For example, the de Rham complex on Lipschitz domains in R.” 
satisfies this property [TB] . 

The abstract Hodge Laplacian is the unbounded operator Lq : kF^ —)> kF^ defined 
hy Lq = d^“^d^+d^_,_^d^ with the domain I?(L§) = {u £ F^DF^* | du S V^_^_^,d*u S 
F^“^}. In the following, we drop the index k when it is clear from the context. For 
example, Lq = dd* + d*d. It is known that for any / £ W^, there exists a unique 
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u G D{Lo) such that 

Lqu = / mod u -L ij. 

Let Kq : / I—> u be the solution operator above. It is known that Kq is self-adjoint 
and compact as a map W ^ W. 

The mixed discretization of the abstract Hodge Laplacian is well-understood. 
The problem above can be formulated in the mixed weak formulation: given / G W, 
find {a,u,p) G x x .5^ such that 

(ct, r) — {u.dr) =0, Vr G , 

(2.1) {da,v) + {du,dv) + {p,v) = {f,v), 

(u, g)= 0 , 

For each k, let be a sequence of discrete subspaces of indexed by h. Vjf is 
called dense if 

VuGV^, lim inf llu — uIIt/= 0 . 
h^OvGVh 

Clearly, density is necessary for to be good approximations of However, 
it is known that this alone is not sufficient for the Galerkin projection to be a 
convergent method. The key additional properties are the subeomplex property that 
dVjf C and the existence of V-bounded cochain projections, that is, bounded 

projections tt* : —>• Vjf satisfying dir^ = iThd. The main result in FEEC states 

that the Galerkin projection of system using dense discrete subspace admitting 
bounded cochain projections is stable. More precisely. Theorem 3.8 of [2] states that 
the bilinear form associated with the system (EID 

Bo{{a,u,p), {T,v,q)) := {a,T) - {u,dT) + {da,v) -K {du,dv) {p,v) - {u,q), 

satisfies the inf-sup condition on x x .5^)^. This implies quasi-optimal 

convergence rates in the H-norm. Further if the bounded projections can be ex¬ 
tended to W-bounded cochain projections, that is, the extension to tt^ : —>■ V/j 

exists and ||7r((||w->-w are bounded uniformly in h, then we get decoupled error es¬ 
timates estimates for each variable in the VF-norm. To state this precisely, we need 
some notations. We use a <b to express that a < Cb for some generic constant C. 
Following [2], we let 

^0 = 11(1 ~ 7r/t)iFo||w'“-s.w''j ho = ||(.^ ~ 

(2.2) ho = lOaxlll (I — nh)dKo\\iYk-j^\Yic-j+i , || (/ — 7r;j)(i*iFo|| }; 

ao = Vo +^o+ ho, 

where Psj is the VF-orthogonal projection from IT to Jo C IF. 

All these quantities converge to 0 as /i 0 due to the compactness and density 
assumption. For example on smooth or convex polyhedral domains for the de Rham 
complex, it is known that r]o = 0{h), 6o = po = 0 (/i’'+^), where r 

is the largest degree of complete polynomials in Vh- We use the following notation 
for the best approximation: for w G V^, 
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Theorem 3.11 of [2] bounds the Vb-norm error of each variable in terms of best 
approximation errors: 

Ik - CT/ill < E{a) + r]oE{da), 

\\d{a - ah)\\ < E{da), 

(2.3) Ik - Uh\\ < E{u) + {ril + 5Q)[E{da) + E{p)\ + 'qo[E{du) + E{a)], 

\\d{u - Uh)\\ < E{du) + r]o[E{da) + E{p)], 
lb -Ph\\< E(p) + fj,oE{da). 

For example, using 'Pr+iE^~^ x VrE^ and assuming full elliptic regularity, we have 
E{a) = and E{da) = for any r > 1. Therefore, we get the optimal 

rates |k — ah\\ = 0(h’'+^). Similarly, the error rates for all variables can be shown 
to be optimal for all four canonical FEEC element pairs. Let P/j : VF —>■ 14, be 
the lE-orthogonal projection and K^h ■ —>■ be the discrete solution operator 
Koh ■ f Uh- These convergence estimates can be stated using operators (cf. 
Corollary 3.17 of [2]): 

\\Ko- KohPhWw^w 

\\dKQ — dKQhPhWw^w + lk*^o — <EhKQhPh\\w^w k ^o- 

The FEEC approach not only gives optimal error rates, but also captures the impor¬ 
tant structures of the Hodge Laplacian. We have the discrete Hodge decomposition: 

The map tt^ ensures that 3ii C 3 and fB/i C 18, but the discrete spaces and 
18)) are generally not subspaces of their continuous counterparts. At the continuous 
level, we have 

(2.5) V/eW, f = dd*Kof + d*dKof + Psjf 

that is, solving the Hodge Laplacian problem with data / leads to the Hodge 
decomposition of /. At the discrete level, similarly, we have: 

(2.6) V/G14, f = ddlKohf + dldKohf + Pf,f€^h®^l®^h- 

This orthogonality is the key ingredient in deriving decoupled W-norm error esti¬ 
mates for each variable above and plays an important role in our analysis as well. 

3. Main results 

We start by identifying with its dual and form a Gelfand triple fl I 4 * C 
Wk ^ (yk PI extend Lq = dd* + d*d to an operator VnV* ^ (VD V*)' 

which is more suitable for studying perturbations: for all u, u G fl V^, 

{Lqu,v) := {d*u,d*v) + {du,dv). 

Our main operator L : 0 F)* —>■ (H^ 0 F)*)' is obtained by perturbing each 

abstract differential: for all u,v £ V'^ (1 , 

(3.1) {Lu, v) := ((d* -I- l 2 )u, (d* -I- ll)v) + ((d + li)u, dv) + {Udu, v) + {I 5 U, v), 

where k : W ^ W are bounded linear maps between appropriate levels for i = 
1,..., 5. More succinctly, we write, 

P = (d © /3)(d © k)i^ © d (d © li^u © l 4 d'u © l^u. 
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The grouping is convenient for the mixed form later. We prove in Lemma l4. II of the 
next section that this L satisfies a Carding inequality. Then standard techniques in 
elliptic PDE theory imply that L is invertible up to some arbitrarily small pertur¬ 
bation to Z 5 . Therefore, generically, it is reasonable to assume that L is a bounded 
isomorphism. 

Let D be the natural domain on which L maps —>■ W^: 

(3.2) := {u e n v; I {d* + l 2 )u e {d + h)u e V^T+i}. 

Our main perturbed problem is: given / G IT*, find u G D such that 

(3.3) Lu = /. 

We reformulate it in the mixed form: given / G IT* find (cr, u) G T*“^ x T* 
satisfying 

(3.4a) (a,T) — (u,dT) — (l 2 U,T)= 0 , Vr G T*“^, 

(3.4b) ((d + l 3 )(T,v) + ((d + li)u,dv) + (l 4 du,v) + (l 5 U,v) = (f,v), Vu G T*. 

The first equation (I3.4al) is equivalent to u G T^* and a = (d* + l 2 )u. The second 
equation (I3.4al) is equivalent to {d + li)u G and d*{d + li)u = f — {d -\-— 
lidu — I 3 U. Hence, if (cr, u) solves (I3H1) . then u solves ( 1 ^ . Therefore, it makes 
sense to use this mixed formulation to solve our problem. 

In addition to the assumption that L is a bounded isomorphism, we need more 
regularity assumptions on li to ensure that L“^(1T) C D so that (13.31) has a solution. 
One of our main tool for analyzing the discretization is the duality argument, where 
the dual problem L'z = g has to be solved as well. Here L' : T* 014 * —>■ (T* OI 4 *)' 
is the dual of L. We collection the conditions under which all these continuous 
problems are well-posed in a theorem: 

Theorem 3.1 (Continuous well-posedness). Let (IT, d) be a Hilbert complex with 
the compactness property. Suppose L defined in equation dSH) is a bounded iso¬ 
morphism and 

(3.5) {d* + l2)L-\W) {d* + 1;){L')-HW) 

Then both the perturbed vrohlem \S.S\ and its mixed formulation \3.4\ are well-posed. 

The proof is given in section |4l The regularity assumption is very mild, without 
which the perturbed problem does not even make sense. The solution operator to 
the dual problem L'z = g will be used frequently, so we give it a name. Let 

K = {L')-^ : (T* n V^y ^ T* n Vy. 

Our first discrete result is the following fundamental theorem on mixed methods 
for problems perturbed by lower-order terms. 

Theorem 3.2 (Discrete stability). Under the assumptions of Theorem \3.1[ sup¬ 
pose further that ||(dd* -I- d*d)K\\w^w is bounded and the following operators are 
compact IT —>■ IT; 

dU^K, {lld-l*2 d*)K. 

Let Vy he a sequence of dense subcomplexes of T* admitting W-hounded cochain 
projections. Then the Galerkin projection of the mixed system dSHl) using the pair 
X Vy is stable in the sense that there exist positive constants ho,Co such that 
for any h G (0,/iq], there exists a unique discrete solution {ah,Uh) G x 

satisfying (13. 4p for test functions in x Vjf and that ||cr/i||y -|- ||Mft,||y < CoH/H- 
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The proof is nontrivial and is given in Section [5l For the de Rham complex, 
suppose all the lower-order terms are multiplication by smooth coefficients and the 
domain has (1 -I-e)-regularity that \\Kf\\ifi+e < ||/||i 2 for e > 0. By definition, 
L' = (d + l 2 ){d* + l^) + {d* + ll)d -j- d*l\ + and K = we have 

{dd* + d*d)K = 1- {lid* + dll + Kd + d*ll + 11+ l^QK 

is bounded —>■ L^. The compactness assumptions are satisfied due to the com¬ 
pact embedding of Sobolev space into for s > 0. This proves the statement 
on the stability of mixed discretization of Hodge Laplacian in the introduction. 

It is well-known that stability guarantees optimal error rates in the energy norm: 

Corollary 3.3. Under the assumptions of Theorem \S.2\. if h < ho, the unique 
discrete solution {ah,Uh) satisfies 

Ik - (ThWv + Ik - Uh\\v < ||(.f - 7r/i)cr|k + ||(^ - TTh)u\\v- 

As is common for mixed methods, the energy norm estimate is crude because 
it couples errors of different variables. For example, for the unperturbed Hodge 
Laplacian solved with the FEEC pair 7^r+iA*+^ x PrA^, the convergence rate in |k|| 
is in fact h^ higher than that in |ku|| but is lumped together with it. Our next finer 
discrete result gives the decoupled bF-norm estimates for each variable similar to 
estimates (1^ for the unperturbed problem. For that, we need more assumptions. 
To simplify the bookkeeping, we define some approximation quantities: 

<5 = max{5o, ||(/ - Ph)llKl ||(/ - Pf,){ll - i;il)K\\, ||(/ - Ph)P<sllKl}, 
(3.6) r] = max{7?o, /tq, <5, ||(/ - Ph)dllK\\, ||(/ - Ph){lld - Z^dk^lL 

a = 5 + rf + po- 

where r]o,do,+o are defined in equation (12.21) and all the operator norms are in 
II • llvu-s-w- As, before, due to the compactness assumptions and density, all S,r],a ^ 
0 as —>■ 0. 

Theorem 3.4. In addition to the assumptions of Theorem \S.‘A assume that 

mid-iid*+ ii-iiii)K\\w^w 

is bounded, then we have the following improved error estimates: 

Ik - <^h\\ < E{a) + {ii + Xi5\/a)E{da) 

+ (X24 + X3V + Xby/d)E{u) + ixsa + Xib'/a)E{du), 

\\d{(r - ah)\\ < E{da) + {xoa + X4 + X 5 V)E{du) 

+ (X45 + XsV + X2X3 )E{u) + (X3 -b X4V)E{(t), 

Ik - rt/ill < E{u) + T]E{du) + r]E{a) -b (a -b X 43 y/a)E{d(j), 

Ikk - Uh)\\ < E{du) + r]E{da) + xi345E{u) + {x3 + xi45i?)A(cr), 

where Xi...j denote the presence of lower-order terms. For example, Xi 25 = 1 k 
k 0 or ^2 k 0 or lo 0, and xi 25 = 0 otherwise. 

The proof is subtle and is given in Section |6l 
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Corollary 3.5. Suppose a Hodge Laplacian problem satisfies full 2-regularity: 
II^/IIh »+2 < ||/||h», foralls> 0 , 

(for example, on a smooth domain), then the error estimates for the discretization 
using FEEC elements are given by Tabled 

Froof. Following the discussion after equation (12.21) . we have po = 0{h), 6o = 
Q(^mm( 2 .r-i-i))^ /ig = where r is the largest degree of complete polynomials 

in Vh- Using the best approximation estimates for FEEC elements, we see that p = 
Oih) and <5 = as well. Plugging these and the best approximation 

estimates into Theorem 13.41 we get the rates in Table [T] □ 

4. Well-posedness at the continuous level 

In this section, we establish well-posedness results for the continuous problem 
and its mixed formulation. 

4.1. Well-posedness for the primal form. Eirst, we prove that the perturbed 
bounded operator is almost always an isomorphism. 

Lemma 4.1. Let {W^,d) be a Hilbert complex having the compactness property 
with domains . Let L be defined as in dSTj. Then, L-\- XI has a bounded inverse 
for all X G C except at a discrete subset (so at most countable). 

Froof. Let M = max^ ||li||vi/^w and 7 = 4M^ -|- M -|- 1/2. Then L-GjI is coercive 
on U'=nt 4 *: 

{Lu, u) 7 (u, u) > (l/2)[(u, u) -b {d*u, d*u) -\- {du, du)]. 

The compactness property ensures that / : n I 4 * ^ {V^ n V/*)' is compact, 

which makes I{L -\- 7 /)“^ compact on {V^ D V^)'. Spectral theory then implies 
that I fJ.I{L 7 /)“^ has a bounded inverse for all p G C except at a discrete 
subset. Then composing with the bounded isomorphism L-\-jI on the right proves 
the claim. □ 

In particular, this shows that either L is invertible oi L el is invertible for any 
small enough nonzero e. 

Then, we prove the well-posedness of our main problem. 

Lemma 4.2. Suppose L defined in equation dSlD is a bounded isomorphism and 
{d* -b l 2 )L~^{W) C V^~^. Then L~^{W) C D, where D is defined in (j3.2ll . In 
particular, problem dll) has a unique solution. 

Proof. Since L is already an isomorphism, we only need to show that L~^{W) C D. 
For any / G W, let u = L~^f. Then by definition, 

{{d*+l 2 )u,{d*+ll)v) + {{d+h)u,dv) + iUdu,v) + {ku,v) = (/,u), Vue 

By assumption, {d* -b ^ 2)11 e . Thus, we have, 

((d -b ii)u, dv) = (/ - (d -b h){d* + i 2 )u - hdu - I 5 U, v), Vu e n u;. 

There is no d*v in the above. By the density of n I 4 * in W^, we conclude that 
the above holds for all v G V^. Hence, (d -b li)u e Thus u G D proves the 

claim. □ 
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We then prove a similar result for the dual problem. Let D' be the natural 
domain on which L' maps —>■ 

D' := V: I {d* + ll)u € (d + ll)u S 

Using the same argument, we get, 

Lemma 4.3. Suppose L defined in equation en) is a bounded isomorphism and 
{d* +q){L')-^{W) C Then {L')-^W) C D'. 

4.2. Well-posedness for the mixed form. We then turn to mixed system (|3.4I) . 
Its associated bilinear form B : (V^~^ x —>• R is: 

(4.1) B{{a, u), (r, u)) := (cr, r) + {du, dv) + {da, v) — {u, dr) 

— {hu, t) + {Isa, v) + {hu, dv) + {Udu, v) + {fiu, v). 

We call this bilinear form well-posed if and only for any {g, /) G x I 4 ', there 
exists a unique solution {a,u) G x satisfying: 

B{{a,u),{T,v)) = {g,T)v'xv + {f,v)v'xv, '^{t,v) G x V'^. 

From the discussion after eauation l3.41 we see that the well-posedness of B implies 
that L~^{W) G D. But it also implies the well-posedness of the dual mixed problem: 
given any {g,f) G x V^, find {^,z) G x satisfying 

B{{p,w), {^,z)) = {g,p)v'xv + {f,w)v'xv, '^{p,w) G x V'^. 

A similar argument shows that the well-posedness of B implies {L')~^{W) G D' as 
well. We collect these results in a lemma: 

Lemma 4.4. Suppose B defined in 63 ]) is well-posed and let L be defined as 
in 631 ). Then L is a bounded isomorphism, L ^{W) C D, and {L') ^(VF) C D'. 

Moreover, the converse is also true. 

Lemma 4.5. Suppose L defined in 631 ) is a bounded isomorphism. Then B defined 
in 631 ) is well-posed if and only if condition ( 1 ^ holds. 

Proof. The only if part is clear from the previous lemmas. We only need to show 
the if part. 

First, we show that B satisfies a Garding-like inequality: there exist positive 
constants a, b, c depending only on such that 

(4.2) B{{a,u),{a,u + ada)) > b{\\a\\v + \\u\\v)^-c\\uf, V(a,«) G x 

Direct computation using Cauchy-Schwarz inequality shows that there exist con¬ 
stants ci,C 2 depending only on HZiHw-s-vu such that 

B{{a, u), {a,u)) > (l/2)(||af + ||d«f) - Ci||uf, 

B{{a, u), (0, da)) > {l/2)\\daf - c^dkf + \\duf + ||uf). 

Multiplying the second inequality by any positive a < l/( 2 c 2 ) and adding it to the 
first inequality, we get the claim. 

Second, fix any (cr, u) G x . We solve a dual problem using u as data: 

let z = cKu and ^ = —{d* -\- ls)z. By assumption, ^ G Direct computation 

shows that 


B{{p,w),{^,z)) = {cu,w), 'i{p,w)€V^ ^xV'". 


FEEC WITH LOWER-ORDER TERMS 


11 


Finally, we add (^, z) to our choice of test functions in the first step and get 
B{{a,u), {a + ^,u + ada + z)) > b{\\a\\v + ||u||v)^- 
Further, from the definition of (^, z), we have 

||cr + ^||v + ||u + Oder + z||y < M||cr||y + HuHy, 

where the constant M depends only on a, ||d(d* + 13 )K\\\y^w, and ||(iitr||iy_s.u/. 
Thus B satisfies the inf-sup condition. Similarly, for any fixed nontrivial (r, v) £ 
X we let u = L~^v and a = {d* + l 2 )u. By assumption, cr € . Direct 

computation shows that 

B{{a,u), {t,v)) = iv,v) > 0. 

It is well-known that these two conditions imply that B is well-posed [3]. □ 

Given these lemmas. Theorem 13.II is clearly true. 

5. Discrete Stability through new projections 

In this section, we prove Theorem 13.21 The idea of the proof is similar to the 
that of the if part of Lemma 14.51 First, due to the subcomplex property, the 
estimate 621) still holds at the discrete level with the same constants a, b, c. Second, 
fix any (cr, u) G x Vu- We can again solve a dual problem using u as data: 

let z = cKu and ^ = —{d* + l^)z. Then we have ^ G ||^||y -I- \\z\\v ^ ||j*|| 

independent of h, and 

B{{p, w), (^, z)) = (cu, w), V(p, w) G X F''. 

But we can no longer add (^, z) to our choice of test functions because (^, z) is not 
discrete. In the rest of this section, we construct a discrete pair {^h, Zh) G F^“^ x F^ 
such that 

UhWv + \\zh\\v < IICIIu + ll^llu uniformly in h, and 

^ ^ |i?((p,u;),(e-a,^-^.))l<ei.(l|p||u + lkllu)(IICIlu + lklk), 

for all (p, w) G X where e;i —>■ 0 as h —>■ 0. Given such a pair, we can add 

it to our choice of test functions: 

B{{a,u), {a + ^h,u + ada + Zh)) > {b - ceh){\\a\\v + ||u||u)^- 
Further, there also exists M > 0 bounded uniformly in h, such that 

lk + 4l|u + \\u + ada + Zh\\v < ^(Ikllu + I|r||u)- 

Choose a sufficiently small hp such that eh^ < bjc. Then for all h < h^, the bilinear 
form B{{a,u), (t, u)) satisfies the inf-sup condition on x with the inf-sup 

constant bounded uniformly below by (6 — ceho)/M. Since F^^“^ x is of finite 
dimension, this establishes the well-posedness. Thus Theorem 13.21 is proved. 

An obvious choice for {^h,Zh) in (15.11) is the elliptic projection given by 

B((p, w), {^H.Zh)) = B{{p, w), (^, z)), V(p, w) G X V^. 

Then eh = 0. But since we have not proved the well-posedness of B on the discrete 
level, we neither know a discrete solution exists nor can we show the uniform 
estimates. The next most obvious choice is obtained using the elliptic projection 
of the unperturbed problem: 

Bq{{p,w,p), {^h,Zh,qh)) = Bq{{p,w), i^z)), \/{p,w,p) G x F^f x 
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Then we have the existence and uniform bounds. But the second estimate in ()5.ip 
fails. In what follows, we develop two new projection operators to correct the elliptic 
projection for the unperturbed problem so that both conditions in (15.111 holds. In 
fact, we do a lot more. Our elaborately chosen {^h,Zh) will not only satisfy (15.IL 
but also have explicit and optimal error rates in quantities like ||^ — \\z — ZhW-, 

and \B{{p,w), — ^hiZ — Zh)\- This is made precise in Theorem 15.31 This result 

also form the basis of the improved error estimates later. Moreover, the two new 
projection operators enjoy many properties making them interesting in their own 
right. 

5.1. Generalized Canonical Projection. Thus far, we have three projections in 
FEEC: the orthogonal projection P^, the cochain projection iTh which commutes 
with d but has no orthogonality property, and the elliptic projection KohPhLo 
which misses the harmonic part and is only well-defined on the subspace D{Lq). 
Here, we introduce a new projection operator H/j : —>■ Vjf given by 

H/t := P3fe -I- dl^KohPhd. 

In the above and for the rest of this paper, for any subspace X of W, we use the no¬ 
tation Px : IT —>■ X for the IT-orthogonal projection. Among other properties, this 
lih satisfies a commutative property generalizing that of the canonical projection 
for classical elements like Raviart-Thomas. 

Theorem 5.1. Suppose (IT, d) is a Hilbert complex satisfying the compactness 
property and are dense discrete subcomplexes admitting W-bounded cochain 
projections. Then H/j is a projection uniformly bounded in the V-norm. Further 
dHh = P< 3 hd. Let rio,ao be defined as in equation (IQ) . Then, for any w € V, 

l|n?iw|| < \\w\\ -\-po\\dw\\, ||(/-n?i)u;|| < ||(J- 7 rft)n;|| -\-r]o\\dw\\. 

Moreover, it satisfies “partial orthogonality”: for any w,v €V, 

|((/- n?,)w,n)| < (||(/- 7 r,i)w|| -l 77 o||dw||)(||(/- 7 r,.)n|| -177o||dn||)-I ao||d^^||||dw||. 

Proof. The stability of the unperturbed discrete problem (12.1|) implies that H/j is 
uniformly bounded in the T-norm. By subcomplex property, we have Phdv = dv 
for V €Vh. Thus, 

UhV = -I- dlKohdv = P:^,,v P^*^v = v, 

showing that H/i is a projection. We know from equation (12. 6 p that for the un¬ 
perturbed problem dd^ATo/i = P'Sh- This proves that dlih = ddl^KohPhd = P<Shd- 
We then prove the first two estimates. Fix any w G V. We split w — lihW = 
(P 3 — P^,,)w {P<s*w — PbjH/iw). The second term can be bounded using the 
error estimates (j2.4p for Koh- 

||P®.w - P<siTlhw\\ = \\d*Kodw - dlKohPhdw\\ < 77o||dw||. 

We then deal with the first term. The subcomplex property ensures f)h (G and 
the cochain property of nh ensures TVhf) C 3ii- These two lead to 

||(P3-P3,HI<||(/-^i.)P3«^II- 

But (J — '!Th)P^w = {I — TThjw — {I — Trh)P<s*w and P<s* = d*KQd. Thus, 

IK-Ps -P5h.)w\\ ^ + Wil - Trh)d*Kodw\\ < ||(/- 7 r,,)w|| -kr?o||dw||. 
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Combining the estimates for the two parts, we get, 

||(/- nft,)w;|| = II (P 3 - P 3 JW + {P<3-w - P< 2 l^hw)\\ < 11(7- 7 r,i)w;|| +r]o\\dw\\. 
By triangle inequality, we get 

||n,,w;|| < llwll + 11(7 - n,,)w;|| < ||w;|| + r]o\\dw\\ 
as well. This proves the first two estimates. Finally, for any w,v G V, we have 

{{I -Uh)w,v) = ((P 3 - P^^)W,V) + ((P®. - P<BlIlh)w, {I - TTh)v) 

+ ((Ps* - P<2l^h)w, TThV). 

The first two terms can be bounded as before. The last term is bounded by 
{d*Kodw - dlKohPhdw,Trhv) = {{Ko - KohPh)dw,TThdv) < ao||fiw||||fiu||, 
where the error estimate (|2.4I) is used again. This finishes the proof. □ 

5.2. Modified elliptic projection. We modify the unperturbed elliptic projec¬ 
tion slightly to accommodate the harmonic forms. 

Theorem 5.2. For any z G D{Lo), let Zh = KohPhLoz + Ps^^Ps^z. Then, 

11^ - Zh\\ < OollPo^^ll, \\d{z - Zh)\\ + \\d*z - dlzhW < ?7o||Po^:||, 
\\Ph{d*dz + dd*z) - {dldzh + ddlzh)\\ < ^o||Po-z||- 

Proof. By equation (|2.5p . we have P^±z = KqLqZ. Thus we have the splitting 

z- Zh = {P^rz - P^rZh) + {P^)Z - Ps^^Zh) = (Kq - KohPh)LoZ + (7 - Psj^)Ps)Z. 

The first term has been estimated by (12.4p . For the second term, since Pj^^^P^j = 
P^^^Pf, and 'Kh'5 C '5h, we have 

Wil - Pf>h)Pf>^\\w < \\{I - '^h)Psjz\\w < Wil - Trh)Ps^\\w^w\\z\\ =^o||2||, 

which proves the first estimate. The second estimate follows from (12.4p directly. Fi¬ 
nally for the last estimate, we use the continuous and discrete Hodge decomposition 

(1231) (EH), we get 

{did + ddl)zh = {did + ddl)KohPhLoz = {P<s^ + P<BOPhLoz. 

Moreover, by definition, Ps^L^z = 0. Thus, 

Ph{d*dz + dd*z - dldzh - ddlzh) = PhLqz - (P©,, -I- P<si)PhLoz = Ps^^Lqz. 

The right-hand side is just ||p —Poll for the unperturbed problem with Lqz as data. 
By eh, 

= ||(Pj-n. -Psj)Loz\\ < 0-f poIIPbPo^II < /io||Poz||, 


which proves the last estimate. 


□ 
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5.3. Projection of the dual solution. We are now ready to construct the discrete 
pair {^h,Zh) satisfying the conditions (|5.ip in the proof of the discrete stability 
theorem. In fact, we prove a stronger result where the first variable p is allowed to 
be in V instead of Vh and derive explicit error estimates. 

Theorem 5.3. Under the assumption of Theorem \3.‘A for any g € , let z = Kg, 

f. = -{d* + ll)z, Zh = KohPhLoz + P^,^Ps^z, and fh = -d^zt - UhlsZ. Then, 

\\z - Zh\\ < ao\\g\\, Mz-Zh)\\<go\\g\\, IIC - 411 ^ Jlllsll- 

Further, for any {p,w) G x , we have, 

{i-^h,z-Zh))\ < [v\\p\\+a\\dp\\ + {fj.o + Xi23V + X5a)\\w\\+X4a\\dw\\]\\g\\. 

Proof. Using the regularity assumption that ||Lo.?^||vu->-iu is bounded, the estimate 
for ll^ — Zh\\ and \\d{z — Zh)\\ follows directly from Theorem 15.21 From the same 
theorem, for we have 

lie - 411 < \\d*z - dlznW + 11(1 - Un)i;z\\ < 770 II 5 II + 114 - U>,)i;Kg\\. 

For the second term, using quantities defined in (13.6p . we have: 

||(/ - Iln)i;Kg\\ < ||(J - ^4^3*^511 + ^I4^3*^ffll < ^bll- 

The last estimate is just a direct computation using the error estimates in Theo¬ 
rem [521 quantities defined in (j3.6L and the Cauchy-Schwarz inequality. □ 

For the proof of Theorem [32] at the beginning of this section, we get {fh, zu) by 
applying this theorem to g = cu G Vh U W. We note that (^, z) here is the same 
as the one defined there. We check that condition (EB is satisfied. First, 

\\zh\\v = \\KohPhLoz + Pf,,,Ps^z\\v < \\z\\v, 

liaik = II - diK.hPhLoz - UhI;z\\v < lld^ll + I4l|v < 1414 +1414, 

where the constants depend only on the stability constant of the continuous and 
discrete unperturbed problem, which is either independent of h or bounded uni¬ 
formly in h. Second, as mentioned before, the compactness assumptions and den¬ 
sity, a, ry, /xq —t 0 as 0. Thus condition (15.111 is verified. This finishes the proof 
of Theorem 13.21 


6. Proof of Improved Error Estimates 

In this section, we prove Theorem [321 To make the notation more compact, we 
use Cu := u — Uh and = {I — Trh)u. Corresponding quantities for a are similarly 
defined. The Galerkin orthogonality equation reads: 

(6.1a) {e,j,T) - {eu,dT) - (heu^r) = 0, yTGV^~^, 

(6.1b) {{d + l 3 )ea,v) + {{d + h)eu,dv) + {l 4 deu,v) + {l 5 eu,v) = 0, VvGVh- 

6.1. Preliuiiuary estiuiates for ||deo-|| aud ||deu||. Optimal estimates for these 
two terms can be obtained directly from the error equations (lO) with carefully 
chosen test functions. 

Lemma 6.1. For any {cr,u) solving (m and {ah,Uh) solving its Galerkin projec¬ 
tion, 

\\de^\\ ^ \\dEa\\ -bxslle^ll -bX4||de„|| -bxslle^ll- 
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Proof. Restricting the test function space to in equation (IB.lbll leads to: 

( 6 . 2 ) Rtg^(de(j l^Ofj ^ 56 ^) 0. 

Thus, dCfj (d P^yf)dCty T Pt^f^dc^^ (d P^yf)dCty T 

Because -k^ maps to IB to IB/j, we have ||(/ — P<Sh)de„\\ < ||(/ — TThjde^jW = ||d£lcr|| 
proving the claim. □ 

Lemma 6.2. For any (cr,u) solving dSl and {ah,Uh) solving its Galerkin projec¬ 
tion, 

IMe„|l < \\dEu\\ + v\\de^\\ + Xi45l|e„|| + XsW^aW- 
Proof. Let Vh = P'si^{'KhU — Uh) in the second error equation (IB.lbl) . We have 

(6.3) {dcu, dvh) = -{dca, Vh) - [{liCu, dvh) + [hca + hcu, Vh)] - {hdcu, Vh). 

By discrete Poincare inequality ||n?i|| < ||diiii|| for Vh G IB^, the second term in (16.31) 
becomes 

\{lieu,dvh) + {hea + heu,Vh)\ < (xi5||e„|| + X 3 ||eo-||)||<in?i||. 

Because deu = d{u — tthu + tthu — Uh) = dEu + dvh, the last term in (16.311 satisfies 

lihdcujVh)] = \{l4dEu,Vh) + {kdvh,Vh)\ < {\\dEu\\ + ||i;/i||)||dnft,||. 

We finally estimate the first term in (16.31) . Let v = P<s*Vh. Then d{'KhV — Vh) = 0 
implies iThV - Vh € 3/i, so {vh - v) E {nhV - Vh). Thus, 

\\v - Vh\\ < ||(/ - TTh)v\\ = 11(7 - 'Kh)d*Kdvh\\ < ri\\dvh\\. 

This implies, 

\{dec,Vh)\ = \{decr,v - Vh)\ < r^Hdec.lHIdr'/ill. 

Combining all these estimates and \\deu\\ < IldTJuH + ||c?n/i|| gives the estimate in 
the claim. □ 

6.2. Duality lemma. The optimal W-norm estimates for ||e„|| and ||ecr|l require 
more work. 

Lemma 6.3. Under the assumptions of Theorem for any g € W, 

\{nheu,g)\ < [?7l|e<T||+a|Me^|| + (/io + Xi23i? + X5a)||e„|| + (Mo?7 + Xi2345a)IMe„||]||g||. 

Proof. Given g, let (^, z) and {^h,Zh) be dehned as in Theorem 15.31 We have, 

{^h^u T g) — 7?((^(T; B/iOii), (^, z)) — R( (e^r , 11/^6^), , Zh ) ) EifCcr , Il^e^t) , (c^ , Cz)). 

Galerkin orthogonality dSH) states that B{{ea, e„), (r, v)) = 0 for all discrete (r, v). 
Thus, 

(Jd-hCu, g) — R((0, (7 Bh)eu) , (^ii , Zh)) Eifccr , n^e^t), (e^, Cz)). 

The second term above can be estimated by the last inequality in Theorem 15.31 and 
l|nhe„|| < ||e„|| +r]\\deu\\, ||dn?,e„|| < ||de^||. 

The result is the following bound: 

\B{{ecr,Ilheu), (e{,e^))| < [? 7 ||e^|| + a\\dea\\ + {p-o + X 1231 ? + X5a)||e«|| 

+ {Tod + Xia)\\deu\W\\g\\. 
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We bound the first term before by splitting it into three parts: 

|B(( 0 , (/ - Uh)eu), {ih,Zh))\ < |((/ - Ilh)eu,d^h) + {d{I - Uh)eu,dzh)\ 

+ \{h{I — n?i)e„, ^h) + {hid — dlh)s-m dzh) + {lb{I ~ dlh)eu, Zh)\ 

+ |(^ 4 c?(/ — n?t)e„, Zh)\ ='■ Qi + Q 2 + Qs- 

We have Qi = 0 because and dUh = P<Sh^- The second term is: 

Q 2 < |((/ - n^)e„, lU + lldz + i;z)\ + |((/ - Uh)eu,l* 2 e^ + llde, + 

We note that + Z*(iz + l^z = [—l^d* + + (Z 5 — The first term above 

can be bounded by the regularity assumptions and the last estimate in Theorem l5.ll 
The second term can be bounded using estimates in Theorem 15.31 The result is: 

Q 2 < [{xi 2 V + X 2 S + X5a)||e„|| + Xi25a|Me„||]||5||. 

For Qsi we have, 

Q 3 = \iil - P<Sh)deu,llzh)\ < \{{I - P<3H)deuJlz)\ + |((/- P<b Jde„, 

For the first term ((/ — P<Q^)deu,l 4 z) = ((ie„, (P^ “ P^h)^ 4 z) can be bounded 
using 13.61 The second term can be bounded by |((ie„, ^ 462 )! and Theorem 15.31 The 
final result is, 

Qs < X 4 a|Me„||||g||. 

Combining the all estimates together, we get the estimate in the claim. □ 

6.3. Preliminary estimate for ||e„||. 

Lemma 6.4. Under the assumptions of Theorem \3.4[ we have 

l|e«|| < \\Eu\\ + r]\\ea\\ + Q!||dea|| + ? 7 ||de„||. 

Proof. Let g = in Lemma 16.31 we get an estimate for (Jlh^mdlhGu)■ Then, 

||e„|| < ||n,je„|| + ||(/ - n,j)e„|| < ||n,je„|| + ||P„|| + r]\\deu\\. 

For sufficiently small /i, we hide the ||eii|| term on the right in the left-hand side. 
The result is the estimate in the claim. □ 

6.4. Preliminary estimate for |lecr||. 

Lemma 6.5. Under the assumptions of Theorem 

llcall < ||£^<T|| + (i7 + X45ya)||rfea||-b(x24 + X3?y + X5v^)||e„||-b(x3a + X45Va)IMeu|l. 
Proof. Equation (I6.1al) implies: 

(C(7, Cfj) (cct; {d Tlfi^e^j.') (cct, Lt/^e^j) 

(CfT, (f (/ 2 CU, dlfiOfy') {e-ui dTlfiCfj'). 

The first term above is bounded by: 

|(e.,(/-n^)e.)|< ||e.||||(/-n^H <(||P.||+r?||dP.||)||e.||. 

The second term is bounded using ||n/jeo-|| < ||ecr|| + 1711^6^11 from theorem l5.11 
|(^ 2 e„,n,jea)| < X 2 ||e«||(||ea|| + r]\\dea\\) 

The last term is estimated by the duality lemma. Because P<s^Ilh = P»h, we have, 
dTlfie^j') — {e-wj P'^i^dou^ — Pgj^^de^j) — ( 11 /^ 6 ^, Pg^^^de^r). 
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We apply lOl with g = P<s,^dea and get 

|(n?,e„, 5 )| < [r?||e^|| + aUde^H + (/io + Xi23?7 + X5a)||e„|| 

+ + Xi2S,iba)\\deu\\]\\P<Shdeo\\- 


From equation (16.2L we have, 

WP^ShdeA < Xalle^ll + X4|Me„|| + Xslle^ll- 

Combining all these estimates, we get the estimate in the claim. □ 

6.5. Proof of Improved Error Estimates Theorem. The estimates in Theo¬ 
rem 13.41 are derived from the four preliminary estimates Lemma 16.11 Lemma 16.21 
Lemma 16.41 and Lemma 16.51 Using 1 for quantities which are bounded and e for 
quantities which goes to zero as /i —>■ 0, the four preliminary estimates has the 
following structure: 


\\deA 


0 

1 

1 

1 ' 




'WdEA 

de^ll 

lle.ll 

< 

e 

e 

0 

e 

1 

0 

1 

1 


\\deu\\ 

l|e.|| 

-(- 

WdEuW 

\\EA 

LlleullJ 


e 

e 

e 

0 


U|e.||J 


\\Eu\\ 


For example, we can substitute the first line into the second and hide the e||dei,|| 
term in the left-hand side of the second line assuming h is sufficiently small. This, 
in effect, switches \\deA in the second line to UdEo-H. Due the vanishing diangonal 
and epsilon lower-triangle of the matrix above, this procedure can be repeated to 
eliminate all unknown error terms like ||e„|| and switch them with known error terms 
like llEull. After this linear algebra exercise, we get the estimates in Theorem 13.41 


7. Numerical examples 

In this section, we show through numerical examples that the error rates given 
by Theorem 13.41 in Table [T] are in fact achieved and cannot be improved. 

In 3D, there are four cases of the Hodge Laplace problems for differential forms 
of degree 0,1, 2, 3. The 0-form and 3-form cases lead to the scalar Laplace problem 
in the non-mixed form and mixed form respectively. The numerical results in these 
two cases are well-known [UlIZllH] and will not be duplicated here. We focus on 
the 1-form and 2-form case. 

Let fl = [0,1]^ be the unit cube in K^. The 1-form mixed Hodge Laplace problem 
with natural boundary conditions is: given f G L^, find u G D satisfying: 

(grad +l 3 ){— div +l 2 )u + curl(curl -|-/i)m -I- 14 curlu -|- l^u = /, in D, 
u ■ n = 0, (curlu + hu) x n = 0, on dXl. 

We choose the following smooth function as the exact solution: 

sin Tra; cos ttz 
cos TTx sin ny , 
cos Try sin ttz 


u = 
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with the coefficients: 


ll — 

sin Try sin ttz 
sin TTZ 

sin TTZ 

sin TTX sin ttz 

sin Try 
sin TTX 

j h — 

o to' 
_1 


sin Try 

sin TTX 

sin TTX sin Try 


1 



■ 3 ■ 


'1 

2 

-T 


'10 

0 

o' 

h = 

2 

, I4 = 

3 

-3 

-3 

, k = 

0 

10 

0 


-1 


1 

3 

1 


0 

0 

0 


At the discrete level, since a = (—div+/2)M is a O-form where PrA° = V~AP , we 
only have two element pairs: VrAP x Vr-iAP and VrAP x V~AP. 

On the same domain 0 = [0,1]^ as before, the 2-form mixed Hodge Laplace 
problem with natural boundary conditions is: given f G L^, find u G D satisfying: 

(curl +/3)(curl +l2)u — grad(div +li)u + h div u + hu. in H, 
uxn = 0, (div-|-Zi)u = 0, on 


We choose the following smooth function as the exact solution in this case: 


(cos TTx -I- 3) sin iry sin ttz 
sin7rai(cos iry + 2)sin7r2; 
sin TTx sin 7r2/(cos ttz + 2) 


with the coefficients 


sin TTX 


'1 

2 

-1' 

— sin Try 

) h — 

2 

-2 

0 

0 


1 

3 

1 



'1 

0 

- 1 ' 


■ 1 ■ 


'10 

0 

o' 

h = 

0 

-1 

0 

, 14 = 

2 

5 ^5 — 

0 

10 

0 


1 

2 

1 


-1 


0 

0 

0 


At the discrete level, we have all four canonical pairs: 


7^r+lA^ X VrAP, 


X VrA^, VrA^ X rpA^, VpA^ x VpA^. 


In all the numerical experiments, we obtain a quasi-uniform triangulation of size 
m for the unit cube H by first triangulating it uniformly with an m x m x m mesh 
and then perturbing each interior mesh node randomly within 20% of the mesh size 
1/m in all three coordinate directions. 

All four pairs of the canonical FEEC elements in dimension < 3 of all degrees 
are supported by the open source finite element package EEniCS [12], in which all 
our numerical codes are implemented. 

For example, for the unperturbed 1-form problem, with V 2 A^ x ViA^, we get 


m 

\W-<^h\\ 

rate 

\\d{a-ah)\\ 

rate 

||u- Uftll 

rate 

\\d{u-uh)\\ 

rate 

2 

2.766e-01 


2.362e-f00 


2.244e-01 


1.441e+00 


4 

3.940e-02 

2.37 

8.529e-01 

1.24 

6.961e-02 

1.43 

7.434e-01 

0.81 

8 

4.504e-03 

3.20 

2.312e-01 

1.93 

1.859e-02 

1.95 

3.740e-01 

1.01 

16 

5.208e-04 

2.98 

5.716e-02 

1.93 

4.659e-03 

1.91 

1.868e-01 

0.96 


Thus, for example, ||(t — ahWh'^ = 2.766 x 10“^ on a 2 x 2 x 2 mesh. The rates 
are computed between the two successive errors. The optimal rates of 3, 2, 2,1 for 
cr, dcr, u, du are clear. With an U lower-order perturbation, we get: 
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m 

\W-<^h\\ 

rate 

\\d{a-ah)\\ 

rate 

||u-u?,|| 

rate 

\\d{u-uh)\\ 

rate 

2 

4.190e-01 


3.522e-b00 


2.215e-01 


1.434e+00 


4 

1.259e-01 

1.44 

1.583e-b00 

0.96 

6.773e-02 

1.42 

7.335e-01 

0.81 

8 

3.689e-02 

1.54 

7.393e-01 

0.96 

1.874e-02 

1.61 

3.749e-01 

0.84 

16 

9.968e-03 

1.86 

3.595e-01 

1.02 

4.799e-03 

1.93 

1.868e-01 

0.99 


Clearly the convergence rates for a and da in are reduced by 1 as predicted. 

There are too many cases for us to list all the detailed results. We instead 
only summarize the numerical results here. First, the error rates in Table [T] are 
guaranteed for all FEEC elements for the Hodge Laplace problem. Second, for 
each case with a reduction in the error rates predicted in Table (TJ there is at least 
one Hodge Laplace problem with a certain form degree which can only converge at 
that reduced rate. In this sense, the rates in Table [I] is optimal. However, we note 
that the rates in Table [T] do no represent an upper bound for all possible cases. For 
example, when the I5 term is given by multiplication by a smooth scalar coefficient, 
we do not observe a reduction of convergence rates in the L^-error of a. We also 
observed that for 1-forms in 3D, an /2-term given by multiplication by a generic 
smooth coefficient do not degrade the L^-error rate in a. 

The full numerical results along with the python source code used in FEniCS 
can be found at the companion code repository at 

https://bitbucket.org/lzlarryli/feeclotexp. 

We note that due to the randomness involved (random perturbation applied to the 
mesh), the error numbers will not be exactly the same but very close to what we 
have listed here. 
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